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ABSTRACT 

We present results of a search for bright Lyman break galaxies at 1.5 < z < 2.5 in the GOODS-S field 
using a NUV-dropout technique in combination with color-selection. We derived a sample of 73 LBG 
candidates. We compare our selection efficiencies to BM/BX- and BzK methods (techniques solely 
based on ground-based data sets), and find the NUV data to provide greater efficiency for selecting 
star- forming galaxies. We estimate LBG candidate ages, masses, star formation rates, and extinction 
from fitting PEGASE synthesis evolution models. We find about 20% of our LBG candidates are 
comparable to infrared luminous LBGs or sub- millimeter galaxies which are thought to be precursors 
of massive elliptical galaxies today. Overall, we can show that although BM/BX and BzK methods 
do identify star- forming galaxies at z ~ 2, the sample they provide biases against those star- forming 
galaxies which are more massive and contain sizeable red stellar populations. A true Lyman break 
criterion at z ~ 2 is therefore more directly comparable to the populations found at z ~ 3, which does 
contain a red fraction. 

Subject headings: galaxies: general — galaxies: high-redshift — galaxies: evolution — galaxies: star 
formation — infrared: galaxies 
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1. INTRODUCTION 

Since their di s covery in the early 1990s 
(jSteidel fc HamiltonI Il993f ). Lyman break galaxies 
(LBGs) have been used very successfully to study galaxy 
evolution processes in the high redshift Universe. Using 
the dropout technique, mult i- wavelength surveys have 
built up large samples of LBGs at redshift s between 
z ~ 3 — 8, where the Lyman brea k is accessible with 
solely ground-based data sets (e.g. iSteidel fc Hamilton 
[I993I; ISteidel et al."1995', ^ 1991 1200I iLehnert fc B7^^ 
[2003; Bouwens et al. 2004). The studies have shown that 
at 2; ^ 4, a significant fraction of the ionizing UV fiux is 
emitted by low luminosity LBGs ( e .g. Lehnert fc Breme: 



20031: iSawicki fc Thompson! 120061 : iBouwens etHI 



:emerl 
2003; 



Lehnert et al.ll201Q[ ). Analysis of the stellar populations 



showed that the population of LBGs at z ~ 5 consists 
of significantly less massive and younger galaxies than 
comparable samp les at z ~ 3 (e.g. Verma et al. 2007; 
lYabe et^H l2009l ). I Verma et al.l suggested that these 
LBGs are experienci ng their first signi ficant buildup of 
stellar mass (see also IStark et al.ll2009f ). At z ~ 3, the 
LBG population consists of more massive galaxies that 
have formed stars o ver longer periods (e.g. Shaplev et al. 
I2OOII : iReddv et alll2008l) . A considerable fraction of the 
z - 3 LBGs (- 15%) consists of massive > I Q ^^^r^) 
infrared luniinous galaxies (jRigopoulou et aP 120061 : 
iMagdis et aPIMc ?): these ILLBGs are considered one 
possible class of progenitors of today's massive elliptical 
galaxies. 

Less is known about the LBGs at z ^ 2.5. At this red- 
shift, the Lyman break occurs in the UV, which can only 



be observed from space. First studies using deep images 
observed by the GALEX satellite obtained sample s of 
LBGs at z ~ 1 using FUV dropouts (Burgarella et al.l 
1200a 120071 : iHaberzettl et al. 2009). The majority of 
LBGs (^ 75%) at this redshift are disk-dominated, 
while the rest co nsist of interacting/merging galaxies 
(|Burgarella et al.l[2006> ). 

At 2; ^ 2, a number of groups have worked around 
UV restrictions by using selection techniques which solely 
rely on ground based data sets as a proxy method for 
identifying star- forming galaxies. Th e BM/BX- method 
(^Adelberger et al.l 120041 : ISteidel et al.l 12004 ) is based on 
UnGR color selection. The color cuts were determined 
from a sample of spectroscopically- verified 1.9 < z < 
2.7 star- forming galaxies, which were fitted with dust- 
reddened SEDs. The resulting best-fits were redshifted 
to z = 1.5,2.0,2.5 and use d to calculate characte ristic 
colors; for more details see (lAdelberger et al.ll2004f ). 

A seco nd ground-based c olor method is the BzK- 
selection (jPaddi et al.ll2004al ). The choice of filters cov- 
ers the fiat UV part of the spectra and fianks the Balmer 
break at rest-frame 4000 A. This method is sensitive to 
selecting and distinguishing between star-forming and 
passively evolving galaxies at 2; ^ 1.4. U sing the K20 
data from the GOODS fields, iDaddi et al.l identified 25 
star-forming and 7 old galaxies and reported a high se- 
lection efficiency with only 13% contamination by inter- 
lopers. 

Follow-up studies showed that the optically (BX) se- 
lected star-forming galaxies at 2: ~ 2 have typical stel- 
lar masses on the order of 3 x lO^^M© and typical star 
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formation rates of ^ 50 yr~^ (e.g. iReddv fc Steidell 
I2QQ4I : iShaplev et all I2QQ5D . Space densities are - 3 x 
10~^ Mpc~^ with a correlation length o f rp = 4.2 ± 
0.4/i-^ Mpc (Reddy et al. 2005; Adelberger^IeLalJ|2005), 
indicating that typical star-forming galaxies at z ~ 2 are 
hosted by ^ 2 X lO^^M© dark matter halos. 

Studies of t he opti c al- NIK selected star-forming galax- 
ies by D addi et al.l ()2005f ) show that a typical BzK 
galaxy at 2; ^ 2 is likely to be a ULIRG with Lir ~ 
(1 - 2) X lO^^L© and star formation rates SFR - 200 - 
3OOM0 yr~^. The space density of these galaxies is 

(1-2) X 10""^ Mp c~^, which is 2-3 times higher than at 
z ^ .1 and 1 (e.g. ISanders et all 120031 : iLe Floc'h et"all 
l2005f ) and makes them relatively common objects at 
z ^ 2. Their ULIRG nature is also supported by morpho- 
logical studies using HST imaging (Daddi et al. 2004b), 
indicating that many of the 5zi^-selected galaxies show 
hints of merger events driving strong star formation 
events. 

Dy namical studies of star-forming galaxies at z 2 
(iFors ter Schreiber et al. 2009), together with the results 
from Burgarella et al.l (120061) , indicate that between red- 
shifts oi z ^ 2 and 2; ^ 1 is when disk galaxies with 
significant masses start to contribute to the co-moving 
density of galaxies. It is also the p hase w hen the star 
formation intensity reaches its peak (iHopk ins 2004, and 
references therein). 

From these solely ground-based selection techniques 
and the development of new instruments like the IFU 
SINFONI at the VLT and OSIRIS at Keck, we have 
learned much about star formation at the critical epoch 
z ^ 2. However, since the selection processes are differ- 
ent to the techniques used for true LBGs at 2: ^ 3, com- 
parisons between the samples are not straightforward. A 
consistent selection method using the Lyman break, in- 
dependent of intrinsic galaxy properties, should be used 
to select consistent samples of star-forming galaxies at 
z 2. Observations require the space-based facilities 
like GALEX or the UV-channels of the newly installed 
WFC3 on b oard HST . First studies include .Hathi et al. 
(|2010[ ) and iLv et al.l The latter analyzed a very deep 
(^140ksec) GALEX NUV image of the Subaru Deep 
Field fSDF lKashikawa et alll2004f ) and other filters to se- 
lect true LBGs at redshifts 1.6 < z < 2.7, and compiled a 
sample of ~8000 LBG candidates down toV = 25.3 with 
NUV ^ 27, including 24 spectroscopically confirmed at 
1.5 < z < 2.7. For z ^ 2 LBGs, they reported a fac- 
tor ~ 1.8 higher summed luminosity density compared 
to z ^ 2 BX and 2: ~ 3 LBGs, and an increase in the 
luminosity density by a factor of 3-6 between z ^ 5 and 

Z r^ 2. 

iH athi et all examined HST+ WFC3 images (F225W, 
F275W, and F336W) of a 50 arcmin^ field of the 
GOODS-S field, as part of an extended multi-color sur- 
vey. The group used [/225 — , ^^275 — , and [/sse— dropouts 
down to AB = 26.5 combined with color criteria to select 
473 LBG candidates at z ~ 1.7, 2.1 and 2.7 respectively. 
The fitted Schechter function parameters agree well with 
predictions from other observations with respect to red- 
shift. 

These two studies demonstrate the feasibility of select- 
ing samples of true LBGs at 2; ~ 2 consistent with higher 
redshift z ^ 3 LBG samples. Although they are a first 



step in better understanding this crucial population of 
star forming galaxies at an epoch formally known as the 
"redshift desert", little currently is known about their 
properties. The studies either lack data cov erage espe- 
cially in the NIR a nd mid-IR dLv et al.ll2QQ9l ) or field of 
view available (H athi et al.l l2010f). 

In this paper, we present first results from a 
dropout+color selection of a sample of true 2; ~ 2 LBGs 
based on deep GALEX observations of the GOODS-S. 
Our specific purpose is to target a rarer, brighter sam- 
ple of LBGs than found in the small, deep survey of 
iHathi et al.l (pQlOl ). In addition , it co mplements the 
wider, shallower study of lLv et aH (|2009f ) because we use 
a wider range of filters and target a different field, to min- 
imize the effects of cosmic variance. We will show that 
we can obtain a sample of bright z ^ 2 LBGs similar to 
the higher redshift population observed at 2; ^ 3. Here 
we report first results of estimates of ages, stellar masses, 
SFRs, and extinction for our purely UV-selected 2; ~ 2 
LBG sample. Using the deep GALEX UV data available 
for the CDF-S, we were able to utilize the complete size 
of the GOODS-S field which offers a wide range of data 
from the UV to mid-IR and also tripled the Hat hi et al.l 
sample size. 

The paper is organized as follows: §2 summarizes the 
data used, §3 describes our search and defines the se- 
lection criteria, and §4, §5, and §6 present photomet- 
ric redshifts, SED modeling and the results. Our sum- 
mary and conclusions are in §7. Throughout the pa- 
per we use a cosmology with Hq = 71 km s~^ Mpc~^, 
= 0.3,andr^A = 0.7. All distances quoted are co- 
moving. 

2. DATA 

We base our study on GOODS and supporting surveys, 
which offer a wide variety of ground- and space-based 
imaging and spectroscopic data sets, all centered on the 
CDF-S and HDF-N. We concentrated our studies on the 
CDF-S, since this field offers the best data coverage. See 
Table [T] for a complete listing. 

2.1. GALEX data 

We used very deep GALEX FUV and NUV imagery, 
which is part of the GALEX Deep Imaging Survey (DIS). 
The data set consists of a 87 ksec FUV and a 89 ksec 
NUV exposure centered on the CDF-S (Martin e t al.l 
I2OO5I ) covering a ~1 deg^ field of view. Without de- 
convolution, the data (FUV and NUV) provide a 3 a 
detection limit of mAB=24.5 (MAGJSO from SExtrac- 
tor). Confusion noise in such deep images, with the large 
GALEX PSF (~5.4''), contributes to decreasing the de- 
tection sensitivity of the data. However, since we were 
interested in the detection of LBGs at z ~2, where the 
Lyman break is shifted in the NUV-band, it is the non- 
detections and the detection limits that are important. 
Therefore, our results are not seriously affected by confu- 
sion or overlapping sources, and consequently we did not 
deconvolve the NUV images. We adopted an effective de- 
tection limit where the completeness reaches ~ 70-75% 
at mAB(FUV)-26.2 and mAB(NUV)-26.7 within a 
aperture (Figure [1]). We calculated the completeness us- 
ing simulated images created with skymaker v3.3.3 and 
verifying detection probability using SExtractor. Fig- 
ure [1] has been generated using a magnitude bin size of 
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Fig. 1. — Results of the completeness calculations for the CDFS GALEX data. Left: FUV Completeness distribution showing a level of 
70-75% at a limit of ~26.2 AB. Right: NUV completeness distribution, indicating a level of 70-75% for ~26.7 AB. 



0.125. The process was repeated 100 times, with 50000 
stars per square degree between 18 and 28 mag per itera- 
tion. The errors for the completeness are represented by 
the standard deviations for the single bins. The smooth 
decline in completeness, starting around ~23 mag, cor- 
relates with the number density of the simulated objects 
and can be interpreted as result of confusion due to the 
large PSF of the GALEX data. The steep decrease at 
faint magnitudes (marked by the big crosses in Figure [T]) 
indicates where the sample starts to become rapidly in- 
complete. We chose a threshold of 70 — 75% because 
it is where the steepest decline sets in and thus becomes 
very sensitive to photometric uncertainties, i.e. the de- 
tection or non-detection of sources becomes sensitive to 
even small photometric errors. 

To avoid false detections or non-detections (for exam- 
ple due to objects being merged with bright and/or ex- 
tended neighbors) we visually inspected each candidate. 
The relatively small number of objects made this inspec- 
tion straightforward, and ~10% of the objects had multi- 
ple identifications which could imply mergers or clumpy 
extinction but probably not separate but overlapping 
sources. 

2.2. Optical and NIR data 

We used deep public images from the GaBods survey 
for the object search and color selection of LBG candi- 
dates (Hildebrandt et al. 2005). The data set offers 0.25 
deg^ UBVRI pointings centered on the CDF-S, and was 
observed using the ESO/MPG 2.2 m telescope with the 
Wide Field Imager (WFI). The images reach 5cr detec- 
tion limits in UBVRI for a 2" aperture of 26.3, 27.2, 26.8, 
26.8. For more details, see Table [TJ 

Deep NIR data are available as part of the GOODS 
survey observed using ISAAC at the ESO VLT 
(jGiavalisco et al.ll2004b[ ). We used data release vl.O cov- 
ering a 10^x16' area observed in the J- and i^g-filters 
and reaching 3a detection limits of <25.5 and <25.1 re- 



spectively. The images yield an average resolution of 
(0.44) arcsec in J and (0.47) arcsec in Kg. 

Besides the ground based broad-band imaging, we also 
used high resolution space-based imaging covered with 
ACS on the Hubble Space Telescope, in particular data 
from GOO DS (Giavalisco et al. 2004a) as weh as GEMS 
()Rix et al.l l2004). observed in the Sloan z-band. The data 
offer a pixel size of 0.03^' or ~250 pc at z ^2, and provide 
us with basic morphological information for 1.5< z <2.5 
LBG candidates. 

2.3. Spitzer data 

Additional photometric information in the mid-IR was 
provided by deep Spitzer images observed in the 3.6, 
4.5, 5.8, 8.0 and 24/im bands. The data are part of the 
GOODS survey and were observed as a Spitzer Legacy 
program in 2004; we used the data version vO.3 (Dickin- 
son et al., in preparation). The images cover 10' x 16'. 
Due to the observing strategy for the IRAC bands (3.6 
to 8.0 /im), which are comprised of 2 exposures taken 
6 months apart and rotated 180° against each other, 
the central 1/3 ultra-deep parts of the final images were 
exposed twice as long as the outer super-deep regions 
(Table [T]). The MIPS data were only observed during 
the epoch 1 observing campaign and correspond to the 
super-deep IRAC images. The Spitzer data cover the 
rest frame NIR for the selected z ^ 2 LBGs, and probe 
the older low mass part of the stellar populations, thus 
better constraining the SED fits. 

2.4. Spectroscopic data 

Redshift information is available for parts of our 
sample from the GOODS-ESO observations with 
VLT+VIMOS (Popesso et alJ[2Q^. We used redshifts 
from the low resolution VIMOS LR-Blue grism data, 
which have a spectral resolution of '^28 A, with a 
5.7 A/pixel dispersion. The spectra cover 3500-6900 A 
(1167-2300 A at 2; = 2). A second set of redshifts used 
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TABLE 1 

Summary of data used. Unless otherwise noted, Mum is a detection limit. 



Survey 


Filter 


Telescope 


tot. Exp. Time 

s 


FWHM 
arcsec 


mim 

ABmag 


Comment 


GALEX DIS 


FUV 


GALEX 


87042 


4.5 


27.0^ 


Martin et al. 2005 


GALEX DIS 


NUV 


GALEX 


89199 


5.6 


26.5" 


Martin et al. 2005 


GaBoDS 


U 


ESO/MPG 2.2 m+WFI 


78891 


1.01 


26.3^ 


Hildebrandt et al. 2005 




B 




69431 


0.98 


27.2^ 


Hildebrandt et al. 2005 




V 




104603 


0.92 


26.8^ 


Hildebrandt et al. 2005 




R 




87653 


0.79 


26.8^ 


Hildebrandt et al. 2005 




I 




34575 


0.93 


25.0^ 


Hildebrandt et al. 2005 


GOODS/GEMS 


F850LP 


HST+ACS 


18232 


0.03 


26.0 


Giavalisco et al. 2004, 
Kix et al. zUU4 


GOODS 


J 


VLT+ISAAC 


(15500) 


(0.44) 


< 25.5^ 


Giavalisco et al. 2004 




Ks 


data release vl.O 


(20200) 


(0.47) 


< 25.2^ 


Giavalisco et al. 2004 




[3.6] 


Spitzer+IRAC 


82800, 165600^ 


3.6 


-24.0 


Dickinson et al. in prep. 




[4.5] 




82800, 165600^ 


3.6 


-24.0 


Dickinson et al. in prep. 




[5.8] 




82800, 165600^ 


3.6 


-23.0 


Dickinson et al. in prep. 




[8.0] 




82800, 165600^ 


3.6 


-23.2 


Dickinson et al. in prep. 




[24.0] 


Spitzer+MIPS 


82800 


3.6 


-20.5 


Dickinson et al. in prep. 




MOS 


VLT+VIMOS 


(14200) 




UcTio ~ 26 


Popesso et al. 2009 




MOS 


VLT+FORS2 


(39000) 




^850 < 26 


Vanzella et al. 2008 



1.5 cr detection limit used within SExtractor search 
^ 5 cr detection limit within a 2'' aperture 
^ total exposure time for the central overlap region 
^ 10 cr detection limit within a V aperture 



in our study was also available v ia GOODS, and was 
obtained w ith the VLT+F0RS2 (iVanzeha et all I2QQ5L 
I2QQ6L 120081 ). The observations were taken using the 3001 
grism without an order-separating filter, resulting in a 
3.2 A/pixel dispersion and a 9 A resolution at 8000 A, 
spanning 6000-10000 A. For this study we could use a 
total of 1075 redshifts in the field, covering 0< z <5. 

2.5. Photometry 

In the following, unless stated otherwise, we will solely 
concentrate on the analysis of objects found in the 
GOODS-S, the inner part of the CDF-S which is cov- 
ered extensively in all 13 filter bands. Since the smallest 
total area coverage (~140 arcmin^) is offered by the NIR 
data observed by ISAAC on the ESO VLT, we used it 
to determine the area over which we selected our final 
sample. 

Throughout this paper we employ both aperture and 
total magnitudes. We used the aperture photometry for 
photometric redshift estimations and SED fitting. To 
perform aperture photometry we applied the aper task 
from NASA IDL Astronomy User^s Library. We mea- 
sured the fiuxes using apertures with equal radii of 2" in 
the optical and the NIR filter bands, based on double the 
FWHM of 1" in the U-band which functions as our detec- 
tion filter. However, since the GALEX and Spitzer data 
offer much lower spatial resolution, we modified the aper- 
ture sizes accordingly (6'' for GALEX FUV- and NUV- 
bands and A" for the IRAC bands). We did not determine 
aperture magnitudes for the Spitzer 24/im MIPS images. 

To account for aperture losses, we derived aperture cor- 
rections in every filter-band using a set of about 100 stars 
found in the GOODS-S field. The aperture corrections 
and their uncertainties range between 0.98 ± 0.78 for the 
GALEX FUV and 0.14 ± 0.17 for the ISAAC (Ta- 
ble [2]) , and were calculated from the growth curves using 
the mean aperture magnitudes and standard deviations 



TABLE 2 
Aperture corrections 



Filter 


radius 


^^aipercor 


^^^aper cor 




II 


ABmag 


ABmag 


FUV 


3.0 


0.98 


0.78 


NUV 


3.0 


0.81 


0.48 


U 


1.0 


0.35 


0.17 


B 


1.0 


0.28 


0.07 


V 


1.0 


0.27 


0.07 


R 


1.0 


0.23 


0.08 


I 


1.0 


0.31 


0.10 


J 


1.0 


0.14 


0.11 




1.0 


0.14 


0.17 


[3.6] 


2.0 


0.34 


0.10 


[4.5] 


2.0 


0.21 


0.08 


[5.8] 


2.0 


0.17 


0.06 


[8.0] 


2.0 


0.31 


0.17 



of the stellar sample. The total errors for the aperture 
magnitudes were determined by adding the uncertainties 
of the magnitude zero-points and aperture corrections to 
the errors determined by aper in quadrature. For the 
Spitzer IRAC-bands we used the following magnitude 
zero-point errors: Ama.e = 0.05 mag, Am4.5 = 0-05 mag, 
Am5.8 = 0.04 mag, Amg.o = 0.06 mag, which were cal- 
culated from the fiux uncertainties of the fiux zero-points 
given in the Spitzer manual. 

We adopted total magnitudes from the SExtractor 
MAG./ 6'0-magnitudes. These total magnitudes were 
used for the color-selection of the LBG sample (see § [3]) 
and the IR color-magnitude diagrams in § I6.2[ 

3. SEARCH AND COLOR SELECTION 

To select LBG candidates in the redshift range 1.5< 
z <2.5, we first detected possible candidates performing 
SExtractor searches in every single filter band. Specifi- 
cally, we did not extract sources from one image based on 
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Fig. 2. — Lower panels: Selection regions defined by a series of model colors of star-forming galaxies in the two color plane, with models 
indicated in the key. Blue dash- dotted quadrilaterals: Colour selection criteria used for LBG candidates at 1.5 < z < 2.5 (the largest 
selection regions, which include all other selection regions). Blue solid quadrilaterals, left panels: BX color criteria (optically selected 
star-forming galaxies for 2.0 < z < 2.6). Blue dashed strip below /hlueward of BX region, left panels: BM color criteria (optically selected 
star-forming galaxies for 1.5 < z < 2.0), BM/BX from Steidel et al. (2004). We select objects for two groups of LBG candidates. (1) Cyan 
diamonds, upper left panel: NUV-dropouts with colors complying with the largest quadrilateral selection box in the B — R vs. U — B 
diagram. (2) Red diamonds, upper two panels: objects at 1.5 < z < 2.0 detected in the NUV-band and possessing colors complying with the 
selection boxes in both color-color diagrams. Blue diamonds: A subsample of 25 NUV-dropouts with spectroscopically confirmed redshifts 
at 1.5 < 2; < 2.5. Black dots: Complete U-band detected sample in the GOODS-S field. 



the detection of an object in another image. Rather, we 
detected objects in every filter band independently, and 
then measured total magnitudes for them. The resulting 
object catalogs were then cross-correlated against each 
other, creating a master catalog including all imaging 
results. For the correlation, we applied matching radii 
according to the spatial resolution of the images used 
(Table [H column 5). As a final step, we matched the 
imaging master catalog with the two spectroscopic cata- 



logs from GOODS (VIMOS and FORS spectroscopy). 

To identify LBG candidates at 1.5 < z < 2.5, we em- 
ployed two different approaches. We limited the size of 
the field searched to the area of the sky covered by the 
ISAAC NIR observations. Since this data set covers the 
smallest area, this step guaranteed to restrict our sample 
only to candidates observed in all 13 filter-bands. Over 
this field we selected all U-band sources brighter than 
25.3 (AB magnitudes) with no detection in the FUV- 
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TABLE 3 

Template SEDs used for deriving color selection criteria 
(Figure EJ. 



model 


star formation law 


age 






[Gyr] 


Irr Bruzual S> 


I, Chariot 2003 constant 


2-3.7 


PEGASE 


burst 


0.5 




exponential deer. 


0.5 




exponential deer. 


1.0 



and NUV-bands (NUV-dropouts i.e. NUV < 26.7 in a 
5" aperture), and compared their locations in a [/ — 5 
vs. B — R color-color diagram to those of model galax- 
ies at 1.5 < z < 2.5 (Figure [2j). The U-band selection 
limit allowed us to detect objects with a Lyman break 
(NUV-U) of at least 1.4 magnitudes. 

We calc ulated mo del galax y color s usin ^ template 
SEDs fromlFioc Ro cca-Volmerangd ([TqQTI . PEGASE) 
and Bruzual fc Ch ariot (2003) , and employed the photo- 
z algorithm HyperZ ( Bolzon ella et al.l l2000l ) to derive 
the model color tracks (Table [3]) for a variety of star for- 
mation laws, ages and dust properties (Figure [2] bottom 
row). The template SEDs were e xtinc tion- adjust ed via 
HyperZ using the reddening law of lCalzetti et all (j2000f ) 
and Ay ranging between and 1.2. We furthermore ac- 
counted for absorption by the Lyg f orest on the SEDs 
using the method from Madau) ()1995l ) within HyperZ. 

The location of the model galaxies in the U — B vs. B — 
R diagram defined the color selection criteria (largest 
blue dash-dotted quadrilaterals in Figure [2j) for z ~ 1.5 — 
2.5 galaxies which, with the NUV Lyman break, we used 
to select our LEG candidates (Figure [2] top row). The 
color cuts are: 

NUV > 26.70 

B-R> -0.53 

U - B>0.70{B - R) -OM 

B-R<0.31{U-B)^0.56 

U-B< (B-R)^ 1.22 (1) 

We applied a second set of selection criteria for LBGs 
at the lower redshift end (1.5 ^ z < 2.0), for which 
the redshifted Lyman break falls within the NUV band 
(2280 A < A < 2736 A). Those color cuts are (Figure [21 
left panels): 

FUV > 26.2 
NUV-U> 1.50 
U-B > -0.41 

U-B <2.59 (2) 



The NUV-selected LEG candidates had to meet both sets 
of color selection criteria, and be detected in the NUV 
but not in the FUV. The U-B limits in the NUV-U\s. 
U — B diagram were the same as in the U — B vs. B — R 
diagram, while the NUV — U > 1.5 limit represents the 
blue end of the model SED tracks. 

Finally, we cleaned our selected sample for stellar 
contaminat ion by c omparing it with a stellar catalog 
([Groenewe gen et aTl [2002). which was created using a 
magnitude-limited, SED fitting approach. This resulted 
in a sample of ~400 stellar candidates in the central 
0.1 deg2 of the CDF-S. About 10% of the NUV-dropout 
and 47% of the NUV-selected galaxies were rejected as 
stars. As a further check, we visually inspected the NUV 
images of the LEG candidates and rejected objects which 
were false negatives (non-detections in the NUV), e.g. 
objects merged with bright neighbors. In total, the clean- 
ing for stellar contamination and the visual inspection for 
false NUV-dropouts reduced our sample size from 171 to 
127 candidates. 

For our subsequent analysis, we concentrate on the 
bright population of LEGs in our target redshift range. 
Therefore, we only included objects which were observed 
in all filter bands between U and 4.5/im. That resulted in 
a sample with uniform sampling across a wide wavelength 
range for SED and photometric redshift determinations. 
This final sample contains 53 NUV-dropout LEG candi- 
dates (see Figure [3] for 2 example images in the available 
filter bands) and 20 NUV-selected LEG candidates down 
toU < 25.3. 

4. PHOTOMETRIC REDSHIFTS 

Using our 13-band photometry from UV to IR 
(GALEX FUV to Spitzer 8/im, but excluding the ACS 
z-band), we estimated photometr ic redshifts with Hy- 
perZ vl.3 (jBolzonella et al.[[2QQQ[ ). The redshift deter- 
mination was done by cross-correlating a set of tem- 
plate spectra to the colors of our sample galaxies. In 
the current version of HyperZ, we us ed a set of five 
template spectra from iBruzual fc Char iot (2003), con- 
sisting of an elliptical, Sc, Sd, Irr, and starburst galaxy. 
In cases where the objects were not detected in a filter 
band, except for the GALEX FUV- and NUV-bands, we 
omitted the non-detection for the redshift determination. 
To account for the Lyman break, the fiux levels in the 
GALEX FUV and NUV band non-detections were set to 
zero during the template correlation. HyperZ can ac- 
count for the effect of the Lya/Ly/3 forest, therefore we 
ap plied the cor responding Lya forest opacity estimates 
of Madau (1995). We accounted for internal reddening 
of the LEG candidates by employing the reddening law 
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of iCalzetti et alJ (|2QQQf ), with Ay set to range between 
and 1.2 mag in steps of 0.2. The resulting redshift dis- 
tributions for galaxy samples selected with the BM/BX, 
BzK and our dropout technique are shown in Figure [H 
The peak at low redshift z ~ in the dropout and NUV 
selected sample (left panel) is mainly due to residual stel- 
lar contamination in the NUV selected sample. Five stars 
were identified when checking the HST images by eye. 
One of the low redshift galaxies, selected with the NUV 
selection method, clearly shows fiux in the FUV filter- 
band. The flux is weak and merged with emission from a 
nearby bright star, leading to a false non-detection. Both 
the star and the galaxy were removed from the sample 
and not considered for further analysis. The second low 
redshift galaxy, selected with the NUV dropout method, 
does not show obvious emission in the UV filter-bands. 
We therefore keep this object in our sample at this stage. 
Our subsequent analysis will use redshift limits to restrict 
the sample further, minimizing the effects on the results 
by foreground or background objects. 

While our NUV dropout selection and the BM/BX 
method allowed for a relatively clean selection of z ^ 2 
star forming galaxies, the contamination with low red- 
shift interlopers is much higher for the BzK selection 
method. 

We calculated photometric redshift uncertainties 
dzoffset = {zphot - Zspec)/{^ + ^spec) (right panel Fig- 
ure [5]). Overall, the rms is Szrms = 0.24, with an aver- 
aged offset of (dZoffset) = 0.13. The fraction of catas- 
trophic outliers {Szrms > 3<j) is 14%. The accuracy of 
the photometric redshift estimation increases with the 
use of our LBG selection method: for our LBG candi- 
date sample, the rms is Szrms = 0.198, with an average 
of {dZoffset) = 0.095, indicating no significant systemic 
offsets. The fraction of catastrophic outliers decreased 
to ~5% (1/20). We therefore conclude that the pho- 
tometric redshift determination enables us to estimate 
unbiased photometric redshift s for our LBG sample. 

5. SED MODELING 

In order to estimate intrinsic properties of the 1.5 < 
z < 2.5 LBG sample, we fitted our 13-band aperture 
photometry to a library of spectral energy distributions 
(SEDs) created by the synthesis evolution model PE- 



GASE (jFioc fc Rocca-Volmerangelll997[ ). We chose the 
PEGASE models mainly because of the possibility to 
model SEDs using consistent chemical evolution which 
constrains the metallicity as a function of time. We think 
this is a more realistic approach tha n using stellar pop- 
ulatio ns with a fixed metallicity (e.g. iBruzual fc Charlotl 
120031 ). Also, due to an otherwise similar modeling ap- 
proach (similar stellar library, similar treatment of the 
TP-AGB phase, etc.) our results are more comparable to 
other studie s of 2; ^ 2 LBGs based onlBruzual fc Charlotl 
models fe.g. IShaplev et al.|[2QQTl 120051 ). 

Our library consisted of several thousand PEGASE 
spectra including star formation histories with constant 
and exponentially decreasing star formation rates. The 
decay times r for the exponentially decreasing SFRs 
range from r = 10 Myr (burst-like scenarios) to r = 
5000 Myr (close to constant sta r formation). 

We modeled the SEDs using a Kroupa et al.' (1993) mi- 
tial mass function (IMF), including consistent chemical 
evolution and following a closed box approach excluding 
gas infall and outflow. Although that it is likely that 
high-redshift star- forming galaxies have outflows (e.g. 
ISteidel et al][20Toh . we decided not to consider outflows 
in the PEGASE SEDs. Again, this allows for a better 
compariso n with previ ously published results based on 
the Bruzu al fc Charlotl models. 

The model flts and the conversion into physical pa- 
rameters (mass, SFR, E{B — V)) were done using an 
IDL based chi-s quare minimization routine. Our rou- 
tine applies the ICalzetti et al.l (|2000 ^ extinction law to 
the PEGASE rnodels for dust correction. We chose 
the ICalzetti et al. law since it has been used by many 
other studies of high redshift g alaxies (e.g. IShap lev et all 
I2001L 120051 : lYabe et al.l 120091 ). We applied extinction 
corrections to the PEGASE model SEDs ranging be- 
tween Ay =0-5 mag in steps of 0.1 mag. The extinction- 
corrected PEGASE models were then redshifted and flt- 
ted following the chi-square minimization approach de- 
scribed bv lGwvn fc Hartwickl (^005). The normalization 
of the measured SEDs was derived including only the 
observed frame optical (UBVR) fllter-bands. 

When possible, we employed spectroscopic redshifts, 
and photometric redshifts otherwise. 

From the model flts we estimated ages, stellar masses. 
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Szrms = 0.198 and a mean offset of (^dzoffset) = 0.095. 



SFRs and dust extinction E{B — V) for LBG candidates. 
Our results only allow for best-fitting SEDs which were 
younger than the age of the Universe. We estimated un- 
certainties for the calculated parameters as the standard 
deviation of the expectation value of each parameter (see 
iNoh et all [20091 for more details). We restrict the final 
sample of LBG candidates to galaxies with photometric 
and spectroscopic redshifts 1.0 < z < 2.5, which includes 
61 objects. 

The SED fitting resulted in a median of the reduced 
values for the complete sample = 1.83 ± 1.70. 
Acceptable fits with a probability of 0.1% (1%) have 
< 3.25(2.5) and make up about 70% (56%) of the 
sample. Larger reduced values appear to result from 
two main sources: 

1) small photometric errors, especially in the NIR and 
mid-IR filter bands, indicating that we are underestimat- 
ing the photometric uncertainties in these bands. 

2) Uncertainties in the photometric redshift estimation. 
We find that reduced x^-values for objects with spec- 
troscopic redshifts (~ 34%) are smaller than for objects 
with only photometric redshifts available. The median 
reduced x^ foi" the objects with spectroscopic redshifts 
is x^ = 1.18 ± 1.53, with about 83% having x^ < 3.25, 
while the median reduced x^ foi" the objects with only 
photometric redshifts is x^ — 2.63 ± 1.70. This indicates 
that uncertainties in the photometric redshift estimation 
resulted in higher uncertainties in the SED fits. Since 
even small redshift uncertainties can have large impli- 
cations for the fit result when including the NUV filter 
(especially for the lower redshift LBGs z ^ 1.5), we de- 
cided to exclude the NUV and FUV photometry from 
the SED fits. 

The majority of best-fitting models are those with 



either constant SFRs (11%) or long decay times r = 
1500 — 5000 Myr (61%) which are in practice nearly con- 
stant. Another 20% of the fitted galaxies are best de- 
scribed by exponentially decreasing SFRs with moderate 
decay times r = 150 — 1000 Myr. Only 8% of the galaxies 
are best represented by exponentially decreasing SFRs 
with short decay times r = 10 — 100 Myr i.e. starbursts. 
All except one of the fitted galaxies are best described by 
dust-reddened SED models. About 84% of the sample is 
best fitted with an extinction correction of Ay > 0.5. 
and 52% with Ay > 1. 

Stellar masses and SFRs are relatively well constrained 
with median uncertainties of about -0.4 dex. We did not 
find an obvious correlation between stellar masses and 
SFRs vs. their x^-values and their errors. However, the 
highest masses (> lO^^M©) seem to be best constrained, 
showing the smallest uncertainties. In contrast, the un- 
certainties for the ages are large, with a median value of 
0.9 dex. Older galaxies in the sample (> 500 Myr) are 
found to have relatively well-constraint ages with a me- 
dian error of ~0.4 dex. Unlike the stellar mass and the 
SFRs, the ages show an increase in uncertainties with in- 
creasing x^- Therefore, most of the uncertainties in the 
best-fitting models lie in the age estimates. 

Examples of successful SED fits are shown in Figure [6l 
The right panel shows an SED for which only photomet- 
ric redshifts were available, while the left panel shows an 
example fit for which the spectroscopic redshifts could 
be used. 

6. RESULTS 

6.1. Dropout selection 

Fro m the GOODS-S spec troscopic data 
(fVanzella et al.l 120081 : iPopesso et all 120091 we were 
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able to obtain spectroscopic redshifts for 25 (20 dropout 
and 5 NUV-selected) of our 67 total LBG candidates. 
About 80% of the dropout selected (16/20) and only 20% 
(1/5) of the NUV selected LBG candidates were found 
to have redshifts in our target range of 1.5 ^ z ^ 2.5. 
However, an additional -15% (3/20) of the NUV- 
dropouts and 60% (3/5) of the NUV-selected objects 
have 1.0 < z < 1.5. Including this poorly explored 
redshift range, inaccessible from the ground for true 
dropout identifications, increased the selection efficiency 
from —80% to 95% for the dropout selected and 20% to 
80% for the NUV selected LBG candidates. 

Comparison with BM/BX. In the last few yea rs, several 
selection methods like the BM/BX- (Steidel et all 120041: 
lAdelberg er et al.ll2004D and BzK-methods (|Daddi et"al] 
l2004ai ib) have been employed, which solely rely on 
ground-based photometry to create samples of star- 
forming galaxies at z — 2. We compared the selection 
efficiency of these methods here to our dropout selection 
which directly probe their Lyman break spectral regions, 
and thus are more directly analogous to the LBG samples 
at z > 3. To avoid effects due to selecting fainter objects 
especially with the BM/BX method, we applied the same 
brightness limits {U < 25.32 AB) as for the dropout 
selection. Figure El shows our Lyman break+color cri- 
teria (large dot-dashed quadrilaterals, left panels) com- 
pared to the BM/BX- method (smaller solid/dashed blue 
quadrilaterals respectively). The overlap is relatively 
small, and consists mainly of model galaxies with rela- 
tively complex star formation histories (e.g. exponen- 
tially decreasing SFRs). The allowed LBG parame- 
ter space increases when including more evolved and/or 
dust-reddened cases - because the dropout method does 
not depend on (optical) colors redward of the rest-frame 
Lyman limit. In summary, only 2/3 of the NUV-dropout 
selected sample of LBG candidates (34 of 53) are de- 
tected with the BM/BX criteria, and only about 60% of 
the BM/BX selected galaxies (26 of 45) would have been 
selected by our color plus NUV-dropout criteria. 

Comparison with BzK. Only 9 LBG candidates (—20% 
of the NUV-dropout selected LBG candidates) satisfy 
the BzK "5 — ZYS.z — Kg^ color requirement. Only 
—25% of the BzK-selected galaxies (11 of 48) were 
found to be NUV-dropouts and about —80% of the 
BzK NUV-dropouts (9 of 11) were located within our 
color-selection region in the U — Bys.B — R diagram 
(Fig.©. 



These results show that a true Lyman break selec- 
tion targets a different population of galaxies, perhaps 
more evolved and/or reddened than the BM/BX-selected 
galaxy population and more similar to the BzK selection. 
Our dropout selection targets mainly bright star forming 
galaxies due to limitations in the sensitivity, especially in 
the K-band. However, compared to the BzK selection, 
the dropout technique reduces the contamination due to 
interlopers and results in a cleaner, more complete sam- 
ple of 1.5 < z < 2.5 LBGs (Figure [71). 

6.2. Masses, SFRs, and Extinction 

Fitting the PEGASE synthesis evolution models, we 
find that our UV-bright LBG sample is best repre- 
sented by ages ranging between 0.004 — 3 Gyr (Fig- 
ure El), stellar masses between 9 x 10^-2 x 10^° M©, 
and SFRs up to 270 M© yr"^ (Figure ^ with a me- 
dian SFR = 50 ± 53M0 yr~^. The sample is relatively 
strongly reddened, with a median Ay = 1.1 ± 0.5. The 
age distribution peaks around 0.1 Gyr (Figure [HI, with 
a possible second smaller peak around 2 — 3 Gyr. We 
s ee a lack of the very young galaxies ^10 Myr found by 
IShaplev et al.l (|2005l ). indicating that our color selection 
compar ed to the BX selected star-forming galaxies at 
z — 2 of lShaplev et al.l is more sensitive to evolved galax- 
ies with moderate star formation rates. However, when 
comparing the stellar mass and SFR distribution, we find 
that the dropout technique is more sensitive to moderate 
masses with a peak aroun d — 10^^ Mc?) and SF Rs higher 
than those at 2: — 3 (e.g. IShaplev et al]l200l[ ). peaking 
at relatively hi^h SFR s of about 100 Mq yr~^. Similar 
tol Verma et al.l (|2007l ). we find that our 2: — 2 LBGs are 
more massive, older, and have similar SFRs compared 
to high red s hift L BGs, for example the z — 5 LBGs of 
lYabe et all (|2009[ ). 

In general, the properties of our dropout-selected LBGs 
based on th e SEP fits are more comparable to LBGs at 
^ - 3 (e.g. 'Shapley et alJ [2001D ^^an are BM/BX and 
BzK-selected galaxies. Our results demonstrate the con- 
sistent selection of the dropout technique over a large 
redshift range, which allows for the compilation of com- 
parable samples to study the evolution of galaxies over 
cosmic time scales. 

6.3. IR Luminous LBGs 
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A small fraction of our LBG candidate sample has 
colors in the R — [3.6] vs. [3.6] color-magnitude dia- 
gram (Figure [TQ|) comparable to infrared luminous LBGs 
(ILLBGs) and sub-mill imeter galaxies at z ~ 3 (e.g. 
iRigopoulou et all l2QQ6f ). The location in the color- 
magnitude diagrams indicates that our LBG sample con- 
sists of about 20% ILLBGs and/or galaxies consistent 
with sub- millimeter galaxies. The estimated median stel- 
lar mass of the IR luminous subsample is 7.9 ± 4.1 x 
10^ M0, and is slightly higher than the median stel- 
lar mass of the total LBG sample (5 ± 4.1 x 10^ M©). 
With a median SFR of 84 ± 62 Mq yr-\ the IR lu- 
minous subsample consists mainly of LBGs with high 
SFRs (>100 Moyr-^), while the distribution for the 



general LBG sample show a tail of LBGs with lower 
SFRs (^30 M0yr~^). The higher median extinction of 
Av = 1.4 ± 0.8 for the ILLBGs {Ay = 1.1 ± 0.5 for the 
complete sample) indicates that light can be subjected 
to opacities well beyond r = 1, even in the optical. The 
implication is that there is still some part of the star 
formation that could be hidden by dust, resulting in an 
underestimate of the SFRs. The high dust reddening of 
the ILLBGs is also supported by the fact that all IR lu- 
minous candidates are found with relatively red colors in 
the color-color diagrams, not falling in the overlap region 
between BM/BX and dropouts {{U - B) = 0.72 ± 0.21 
and {B — R) = 0.53 ±0.27 consistent with r ecent results 
for a small far-IR detected sample by iBurgarella et al.l 
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[2QTTI) . Although our results indicate substantial redden- 
ing due to larger amounts of dust, we have to be care- 
ful when drawing conclusions concerning extinction from 
UV and optical data alone. The age-extinction degener- 
acy and the lack of FIR data, where the thermal emission 
of dust defines the spectral properties, makes any inter- 
pretation difficult at best. Using the results from the best 
fitting SEDs of our sample, we determine a median error 
for the extinction of A Ay ~ 0.5. In all cases we find 
no acceptable fits with Ay >2.3 mag, with the largest 
extinction corrections being found for the ILLBGs. 

We also ca lculated the total IR luminosity following the 
approach of iRieke et al.l ()2QQ9f ), estimating Ltir from 
the 24 jiim flux using a linear approximation. The lumi- 
nosities are pven in L^. emplo ying apparent solar magni- 
tudes using Rieke et alj (|2QQ8f ) and'Engelke et al." (2010). 
We thus estimated total IR luminosities for our subsam- 
ple of ILLBGs of 6 x 10^° - 1 x lO^^L^, with a median 
total IR luminosity of Lttr = 4 x l(r^Lc?) . Using the 
conversion reported by (jNordon et al.l 120101 ) this trans- 
lates into a factor of two lower median SFR of about 
40 Mq yr~^ for the ILLBG sample than derived from 
the SED model fits. 

Compared to LBG samples at higher and lower red- 
shifts with similar mass and UV luminosity ranges, we 
find that our sample consists of galaxies which are mod- 
erately higher extinction. Studies at z ~ 5 and z ^ 3 
estimate median extinctions of Ay ^ 0.3 — 0.9 (e.g. 
Verma et al. 2007; Yabe et al. 2009) and Ay = 0.5 (e.g. 
Shaplev e t al. 2001) respecti vely. At z - 3 the ILLBG 
sample of Rig opoulou et aH (|2006i) shows a median ex- 
tinction of Ay r^l.O. Although slightly higher than the 
general samples of LBGs at this redshift, their reddening 
is lower than the median extinction of Ay = 1.4 found in 
this study. Extinction values {Ay = 0.8 — 1.8) similar to 



ours have been estin iated for a 8iJ,m an d 24/im detected 
sample of LBGs bv iShim et al.l ()2007[ ). However, their 
sample consists mainly of ULIRG-like, massive LBGs 
with stellar masses of up lO^^M© and IR luminosities 
of Ljji > IO^^Lq. Thus our selection, which in general 
selects z ^ 2 LBGs, nevertheless favors selecting galaxies 
with higher extinction compared to z > 3 LBGs. This 
is because we require them to be detected in all bands, 
including in the IR. 

At lower redshift s, studies revealed LBGs with ex- 
tinction values similar to th e ones found at z > 3. 
At z ^ 1 Nilsson et al.| ()2011l ) estimated extinctions of 
Ay = 0.5 — 2, with a m edian of 0.8 for a samp le of LBGs 
in the GOODS-S field (.Burgarella et al.ll2007n . Similarly, 
at z ~ 0.1 — 0.2, Lyman Break Analogs (L BAs) are found 
to ha ve a median extinction of Ay = 0.8 (Overzier e t al.l 
l2011f ). While overall trends in the extinction with red- 
shift are difficult to determine, it appears that our study 
stands out in selecting galaxies with moderately higher 
extinction on average than other samples at lower and 
higher redshift s. 

Compared with other published examples of IR lumi- 
nous LBGs, we find that our sample consists of galax- 
ies with moderate masses having a median mass of 
^8xl0^Mr:^. This i s comparable to the z ^ 1 LBGs of 
iNilsson et al.l ()2011f ). which have masses between 1.3 x 
lO^M© and 3.6 x lO^^Mr;^. At ^ ^ 3 , ILLBGs appear to 
be more massive. "Rigopoul ou et all (2006), for example, 
found IR bright LBGs with masses between 4 x lO^^M© 
and 4 x 10^ ^M©. Even higher masses of 1 0^^M(7) have 
been reported for 24/im-detected LBGs by iShim etld] 
(j2007[ ). With IR luminosities of Ljr > lO^^L©, these 
galaxes are massive, dusty ULIRGs. The very high mass 
LBGs are relatively rare and might have been missed by 
us due to our relatively small search volume. Also, the 
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Fig. 10. — [3.6] vs R— [3.6] (left panel) color-magnitude diagram using Spitzer measurements for our LBG candidate sample. Black dots: 
All galaxies in our sample with either IRAC 3.6 /im (left panel) or IRAC 8.0 fim and MIPS 24 /xm (right panel) detections. We compare 
our measurements to results for z ~ 3 ILLBGs reported by R igopoulou et al.l (|2000, open diamonds). We also make a similar comparison in 
the [8.0] — [24] vs [8.0] color-magnitude diagram (right panel). Orange diamonds, left panel: LBGs from our sample which are also detected 
in the Spitzer 8jLi.ni- but not in the 24/i,m-band. Open diamonds without open stars, right panel: LBGs from the z ~ 3 LBG sample of 
IRigopoulou et all (|2006l ) detected in 8fim but not 24/im. 



lack of very massive LBGs in our sample might have been 
due to limiting our selection to UV-bright objects, requir- 
ing detections from U to IRAC 4.5/im (observed frame). 
Both of these criteria (UV and IR luminous) favor the 
selection of objects with patchy extinction. Such galax- 
ies possess strong, obscured starbursts. However, the 
extinction is patchy such that light from a UV-emitting 
stellar population, or less obscured star formation, can 
be observed. This can be seen in Fig. [TTl where some of 
the galaxies have a complex morphology. The ILLBGs 
represent the more massive, dusty, heavily star- forming 
tail of the LBG population in our study, but are less mas- 
sive than their high SFR, dusty counterparts at higher 
redshifts. 

These galaxies also might resemble an interesting part 
of the high redshift sources contributing to cosmic IR 
background (GIB) light. At z ^ 2 about 30 -40% of the 
GIB results from LIRGs (Ber ta et alll2010D . which are 
comparable to our ILLBGs. In addition, far-infrared ob- 
servations show that at A > 350/im more than half of the 
contribu tion to the GIB co mes from sources at redshifts 
z > 1.2 (|Devlin et al.ll2009l ). Thus, a significant fraction 
of the SFR at this epoch might be provided by dusty 
IR-bright galaxies. 

6.4. Morphology 

From an inspection of the high resolution HST-AGS 
images by eye, we conclude that the majority of the 
LBG sample consists of galaxies showing indications of 
disk-like structure, while the rest can be considered com- 
pact objects (examples shown in Figure [TT]) . Also, the 
majority of the sample has at least one close (< 3^') 
neighbor or shows disturbed disk structures, such that 
they may be considered merger candidates. For the 
ILLBG subsample, we find that about 80% show indi- 



cation of a disk structure and 20% are found to be com- 
pact objects. About 50% of the ILLBGs can be con- 
sidered merger candidates. Although these results have 
to be taken carefully, since we did not apply detailed 
morphological analysis and a judgment by eye is very 
subjective , they seem to b e in agreement with results 
found by BurgardlaeFaT (2006) for the z - 1 LBGs 
which consist of about 75% disk-dominated galaxies and 
z ^ 1 LIRGs (e g. [Melbourne et al.l l2005l: iBeh et al.l[2005l : 
iHammer et al.ll2007l ). A more detailed analysis of the 
morphology will be presented in a later paper (Haberzettl 
et al. in prep.). 

7. SUMMARY AND GONGLUSION 

We present results of a search for Lyman break galax- 
ies at 1.5 < z < 2.5, an epoch once called the redshift 
desert. Over the l ast few years , studies employing the 
the BM/BX- (Steide l et al]f2003 : 1 Adelb erger et al.l[200l 
and BzK-methods (Da ddi et al." 2004a. b) used observed- 
frame optical color-selection criteria to find and study 
star-forming galaxies in this redshift range - which ex- 
cluded the Lyman break. Such methods identify galax- 
ies which, though star-forming, are not completely con- 
sistent with LBG samples at higher redshift {z > 3). 
The difference is that the dropout method allows for a 
wider range in color redward of the Lyman limit, which 
in practice means including redder, more evolved and 
potentially more massive galaxies. 

Our NUV-dropout plus optical color selection tech- 
nique targeted star- forming galax ies at 1.5 < z < 2.5, 
and is similar to the one used bv iLv et"al] (|2009l ). Us- 
ing the deep publicly available data set of the GOODS- 
S (CDF-S) field covering the UV to mid-IR, we iden- 
tified a sample of 73 LBG candidates. Our selection 
method resulted in a success rate of ~ 80% (~95% for 



LBGs at ^ ~ 2 
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Fig. 11. — Example images showing the morphological classification applied. The image size is 3'' x From left to right, and referring to 
galaxies by catalog number: disk-like (top row): #11800 (dropout-selected LBG), #7312 (dropout-selected ILLBG), #6039 (NUV-selected 
LBG), #7479 (NUV-selected ILLBG); compact (bottom row): #6633 (dropout-selected LBG), #12478 (dropout-selected LBG), #7312 
(dropout-selected ILLBG), #10869 (NUV-selected LBG). 



1.0<z< 2.5). 

The NUV dropout technique ahows for the detection 
of more evolved and/or dust-reddened galaxies not de- 
tected with the BM/BX- method, while also decreasing 
the number of interlopers compared to the BzK-method. 
Although the BM/BX and BzK methods can select much 
fainter objects than we can here, due to the limitation 
of GALEX in sensitivity, they exclude the more evolved 
and/or dust-reddened galaxies completely. The fact that 
40% of the dropout-selected LBGs are not detected by 
the BM/BX and/or BzK methods provides for additional 
UV/optical-selected star formation at z ~ 2, especially 
if hidden by large amounts of dust. 

We find that our z ^ 2 galaxy sample is consistent 
with LBGs at z ^ 3, and therefore allows for consistent 
follow-up studies at redshifts spanning z ^ 1.0 — 7. 

We calculated useful photometric redshifts using Hy- 
perZ and included a small set of spectroscopic redshifts 
for our LBG candidate sample (Figure [5l red filled cir- 
cles). The normalized rms scatter for our LBG candidate 
sample is Szrms = 0.198, with a smaller mean offset of 
(dZoffset) = 0.095. 

Using the combination of spectroscopic (where avail- 
able) and photometric redshifts, we fitted synthetic mod- 
els from a library of PEGASE spectra to our LBG sam- 
ple, and determined ages ranging from 4 Myr to 3 Gyr, 
masses between 9 x IO^Mq to about 2 x lO^^M©, and 
SFRs up to ^ 270 Mq yr"^. 

Our LBG candidates have moderate stel lar masses and 
SFRs consistent with z ^ 3 LBGs (e.g. iShaplev et al.l 
l2001f ). Ages and dust extinction values of our sam- 
ple are also more consistent with LBGs at z 3, in 
comparis on with th e BX-selected galaxies at z ~ 2 of 
(e.g. Sha plev et al.l [2005) . although our sample shows a 
lack of the very young star-forming galaxies found in the 
BM/BX and z ^ 3 samples (Figure [8]). Our sample rep- 
resents the moderate mass range of the UV luminosity 
distribution of LBGs at this redshift. This can lead to se- 
lection biases against objects with much higher masses. 
We therefore should exercise care when comparing our 



results to other surveys. 

We also identified a small fraction of our LBG sam- 
ple (~20%) as potential ILLBGs and/or sub- millimeter 
galaxies using Spitzer 8.0 and 24 jam data. Although 
there is some overlap with the other z ^ 2 LBGs, the 
R- [3.6] vs. [3.6] plot in Figure HO] shows that the ILLBGs 
(blue and orange diamonds) represent the UV-bright end 
of the LBG distribution for a given value of i? — [3.6]. Due 
to their red colors at optical wavelengths (compared to 
the BM/BX selected), most of these galaxies are missed 
by the BM/BX and BzK methods (- 10% detected). 

The ILLBG SED fits result, on average, in higher total 
stellar masses, SFRs, and extinctions than our overall 
sample. 

Compared with higher and lower redshift LBGs from 
the literature, we find that our z ^ 2 sample has an 
average extinction which is higher than LBGs at ^ 3 (e. g 
IShaplev et aPIMnl : iVerma et al.l[2007l : lYabe et al.ll2009l) . 
A similar resu lt holds in c omparison with ILLBGs at 
2: ^ 3 (e.g. Rig opoulou et a l. 2006), with the exception 
of massive, IR luminou s, ULIRG-like L BGs which are 
selected in the infrared (jShim et al.ll200Z) . We find that 
the masses and extinction values are in better agreement 
with results for 2; ^ 1 LBGs 



Nilsson et al.ll201lh . 

The median total IR luminosity of our subsample of 
Ltir = 4x 10^^ Lq indicates that these ILLBGs are com- 
prised mainly of LIRGs. Only - 15% (2 out of 15 ILL- 
BGs) are comparable to ULIRGs with Ltir > lO^^L©. 
From the total IR luminosities, we estimate a median 
SFR half that compared to the one derived from the SED 
model fits. The role of ILLBGs in the overall scheme of 
star formation at z > 1.2 and their contribution to the 
cosmic IR background merits further investigation, as 
ILLBGs could contribute significantly. 

A simple morphological classification shows that the 
majority of our LBG candidates appear disk-like, and 
many show indications for mergers and interactions. 
Similar results were found for the ILLBGs, with 80% 
having disk-like morphologies. T his is consistent with re- 
sults published for z - 1 LBGs (jBurgarella et al.ll2006D 
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and z ^ 1 LIRGs (e g. fMelbQuriie eraIl l2QQ5l: iBell et all 
I2QQ5I : [Hammer et alJl2QQ7l ). 

Overall, our study shows that the selection of star- 
forming galaxies at 1.5 < z < 2.5 using a true Lyman 
break dropout technique allows for the selection of a ro- 
bust sample of LBGs, which can used for comparable 
follow-up studies with LBGs at higher redshifts. A Ly- 
man break selection like we used here allows redder, more 
massive LBGs to be included in comparison to other se- 
lection techniques such as the BM/BX and BzK meth- 
ods. 
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